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Skyrmions recently discovered in chiral magnets are a promising candidate for magnetic storage devices 
because of their topological stability, small size (~ 3 — lOOnm), and ultra-low threshold current density (~ 
10®A/m^) to drive their motion. However, the time-dependent dynamics has hitherto been largely unexplored. 
Here we show, by combining the numerical solution of the Landau-Lifshitz-Gilbert equation and the analysis of a 
generalized Thiele’s equation, that inertial effects are almost completely absent in skyrmion dynamics driven by 
a time-dependent current. In contrast, the response to time-dependent magnetic forces and thermal fluctuations 
depends strongly on frequency and is described by a large effective mass and a (anti-) damping depending 
on the acceleration of the skyrmion. Thermal diffusion is strongly suppressed by the cyclotron motion and is 
proportional to the Gilbert damping coefficient a. This indicates that the skyrmion position is stable, and its 
motion responds to the time-dependent current without delay or retardation even if it is fast. These findings 
demonstrate the advantages of skyrmions as information carriers. 

PACS numbers: 73.43.Cd,72.25.-b,72.80.-r 


I. INTRODUCTION 

Mass is a fundamental quantity of a particle determining its 
mechanical inertia and therefore the speed of response to ex¬ 
ternal forces. Furthermore, it controls the strength of quantum 
and thermal fluctuations. For a fast response one usually needs 
small masses and small friction coefficients which in turn lead 
to large fluctuations and a rapid diffusion. Therefore, usually 
small fluctuations and a quick reaction to external forces are 
not concomitant. However a “particle” is not a trivial object 
in modem physics, it can be a complex of energy and mo¬ 
mentum, embedded in a fluctuating environment. Therefore, 
its dynamics can be different from that of a Newtonian parti¬ 
cle. This is the case in magnets, where such a “particle” can 
be formed by a magnetic texturJi^ A skyrmiorP® is a rep¬ 
resentative example: a swirling spin texture characterized by 
a topological index counting the number of times a sphere is 
wrapped in spin space. This topological index remains un¬ 
changed provided spin configurations vary slowly, i.e., dis¬ 
continuous spin configurations are forbidden on an atomic 
scale due to high energy costs. Therefore, the skyrmion is 
topologically protected and has a long lifetime, in sharp con¬ 
trast to e.g. spin wave excitations which can rapidly decay. 
Skyrmions have attracted recent intensive interest because of 
their nano-metric size and high mobilitj®^. Especially, the 
current densities needed to drive their motion (~ 10®A/m^) 
are ultra small compared to those used to manipulate domain 
walls in ferromagnets (~ 

The motion of the skyrmion in a two dimensional film can 
be described by a modified version of Newton’s equation. For 
sufficiently slowly varying and not too strong forces, a sym¬ 
metry analysis suggests the following form of the equations 
of motion, 

Q X R ofD It -f TTiit -f crF X It = (1) 

Here we assumed translational and rotational invariance of 
the linearized equations of motion. The ‘gyrocoupling’ Q = 


Q e\_ is an effective magnetic held oriented perpendicular to 
the plane, a is the (dimensionless) Gilbert damping of a sin¬ 
gle spin, olD describes the friction of the skyrmion, m its mass 
and R its centre coordinate. F parametrizes a peculiar type of 
damping proportional to the acceleration of the particle. We 
name this term ‘gyrodamping’, since it describes the damping 
of a particle on a cyclotron orbit (an orbit with R cx ^ x R), 
which can be stronger (F parallel to Q) or weaker (antipar¬ 
allel to Q) than that for linear motion. Our main goal will 
be to describe the influence of forces on the skyrmion arising 
from electric currents (Fc), magnetic held gradients (Fg) and 
thermal fluctuations (Fth)- 

By analyzing the motion of a rigid magnetic structure 
M(r,f) = Mo(r — R{t)) for static forces, one can obtain 
analytic fo rmula s for Q,aD,Fc and Fg using the approach 
of ThieldSmai In Ref. an approximate value for the 

mass of a skyrmion was obtained by simulating the motion of 
a skyrmion in a nanodisc and by estimating contributions to 
the mass from internal excitations of the skyrmion. 

For rapidly changing forces, needed for the manipulation of 
skyrmions in spintronic devices, Eq . (p^ is however not suffi¬ 
cient. A generalized version of Eq. (fTTvalid for weak but also 
arbitrarily time-dependent forces can be written as 

= F,M+Fg(w) + Ft;,(a;) (2) 

= Sc(a;)t>s(a;) + Sg(w)VB2(w) + Fth(w) 

Here V{uj) = f e^‘^*R{t)dt is the Eourier transform of the 
velocity of the skyrmion, Vg{uj) is the (spin-) drift velocity 
of the conduction electrons, directly proportional to the cur¬ 
rent, V Bz (uj) describes a magnetic held gradient in frequency 
space. The role of the random thermal forces, Fth(w), is spe¬ 
cial as their dynamics is directly linked via the fluctuation- 
dissipation theorem to the left-hand side of the equation, see 
below. The 2x2 matrix describes the dynam¬ 

ics of the skyrmion; its small-w expansion defines the terms 
written on the left-hand side of Eq. Q- One can expect 
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FIG. 1: When a skyrmion is driven by a time dependent external 
force, it becomes distorted and the spins precess resulting in a de¬ 
layed response and a large effective mass. In contrast, when the 
skyrmion motion is driven by an electric current, the skyrmion ap¬ 
proximately flows with the current with little distortion and preces¬ 
sion. Therefore skyrmions respond quickly to rapid changes of the 
electric current. 


strongly frequency-dependent dynamics for the skyrmion be¬ 
cause the external forces in combination with the motion of 
the skyrmion can induce a precession of the spin and also ex¬ 
cite spinwaves in the surrounding ferromagnet, see Fig. [T] 
We will, however, show that the frequency dependence of the 
right-hand side of the Eq. (2) is at least as important: not only 
the motion of the skyrmion but also the external forces excite 
internal modes. Depending on the frequency range, there is 
an effective screening or antiscreening of the forces described 
by the matrices Sc(a;) and Sg(a;). Especially for the current- 
driven motion, there will be for all frequencies an almost exact 
cancellation of terms from and Sc(a;). As a result 

the skyrmion will follow almost instantaneously any change 
of the current despite its large mass. 

In this paper, we study the dynamics of a driven skyrmion 
by solving numerically the stochastic Landau-Lifshitz-Gilbert 
(LEG) equation. Our strategy will be to determine the param¬ 
eters of Eq. (j^ such that this equation reproduces the results 
of the LEG equation. Section |I^ introduces the model and 
outlines the numerical implementation. Three driving mecha¬ 
nisms are considered: section UII] studies the diffusive motion 
of the skyrmion due to thermal noise, section|^the skyrmion 
motion due to time-dependent magnetic field gradient and sec- 
tion[^the current-driven dynamics. We conclude with a sum¬ 


mary and discussion of the results in Sec. VI 


II. MODEL 

Our study is based on a numerical analysis of the stochastic 
Eandau-Eifshitz-Gilbert (sEEG) equation^ defined by 

^ = 7 Mr X [Beff + bfl(f)] 

- T^Mr X (Mr X [Beff + bfl(f)]) . (3) 

Here 7 is the gyromagnetic moment and a the Gilbert damp¬ 
ing; Beff = — is an effective magnetic field created by 

the surrounding magnetic moments and bfl(f) a fluctuating. 


stochastic field creating random torques on the magnetic mo¬ 
ments to model the effects of thermal fluctuations, see below. 
The Hamiltonian iJ [M] is given by 


H[M] = -J^Mr • (Mr+ae,, + Mr+aeJ 

r 

- A ^ (Mr X Mr+ae,, • + Mr X Mr+ae„ ' Gy) 

r 

- B-^Mr (4) 

r 

We use J = 1, 7 = 1, I Mr I = 1, A = 0.18J for the 
strength of the Dzyaloshinskii-Moriya interaction and B = 
(0, 0, 0.0278J) for all plots giving rise to a skyrmion with a 
radius of about 15 lattice sites, see Appendix [A| Eor this pa¬ 
rameter set, the ground state is ferromagnetic, thus the single 
skyrmion is a topologically protected, metastable excitation. 
Typically we simulate 100 x 100 spins for the analysis of dif¬ 
fusive and current driven motion and 200 x 200 spins for the 
force-driven motion. Eor these parameters lattice effects are 
negligible, see appendix Typical microscopic parameters 
used, are J = 1 meV (this yields Tc ~ 10 K) which we use to 
estimate typical time scales for the skyrmion motion. 

Eollowing Ref. |27] we assume that the held b^(f) is gen¬ 
erated from a Gaussian stochastic process with the following 
statistical properties 

KM = 0 

(KAM'.jis)) = 2a^^Sij6rr'6{t - s) ( 5 ) 


where i and j are cartesian components and (...) denotes 
an average taken over different realizations of the fluctuating 
field. The Gaussian property of the process stems from the in¬ 
teraction of Mj. with a large number of microscopic degrees 
of freedom (central limit theorem) which are also responsi¬ 
ble for the damping described by a, reflecting the fluctuation- 
dissipation theorem. The delta-correlation in time and space 
in Eq. Q expresses that the autocorrelation time and length of 
b^(<) is much shorter than the response time and length scale 
of the magnetic system. 

Eor a numerical implementation of Eq. 0 we follow 
Ref. 122 ] and use Heun’s scheme for the numerical integration 
which converges quadratically to the solution of the general 
system of stochastic differential equations (when interpreted 
in terms of the Stratonovich calculus). 

Eor static driving forces, one can calculate the drift veloc¬ 
ity R following ThielJ^. Starting from the Eandau-Eifshitz 
Gilbert equations, Eq. we project onto the translational 
mod e by m ultiplying Eq. ® with O^Mr and integrating over 
spacdSEH. 


g = 

V = 
-Fc = 

Fn = 


hMo J dr n- {d^n x dyu) 

HMq j dr [dxTi ■ dxTi + dyU ■ dyii)/2 
g X Vs + PVVs, 

MsVB, Ms=Mo Jdr{l-n,) (6) 
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where n is the direction of the magnetization, Mq the lo¬ 
cal spin density, Vg the (spin-) drift velocity of the conduc¬ 
tion electrons proportional to the electric current, and Mg 
is the change of the magnetization induced by a skyrmion 
in a ferromagnetic background. The ’gyrocoupling vector’ 
Q = (0,0, C/)"'" with Q — zbSModTT is given by the winding 
number of the skyrmion, independent of microscopic details. 


III. THERM AT, DIFFUSION 


Random forces arising from thermal fluctuations play a de¬ 
cisive role in controlling the diffusion of particles and there¬ 
fore also the trajectories R{t) of a skyrmion. To obtain R{t) 
and corresponding correlation functions we used numerical 
simulations based on the stochastic Landau-Lifshitz-Gilbert 
equatiorP^. These micromagnetic equations describe the dy¬ 
namics of coupled spins including the effects of damping 
and thermal fluctuations. Initially, a skyrmion spin-texture 
is embedded in a ferromagnetic background. By monitoring 
the change of the magnetization, we track the center of the 
skyrmion R{t), see appendix|^for details. 

Our goal is to use this data to determine the matrix 
and the randomly fluctuating thermal forces, J’th(tu), which 
together fix the equation of motion, Eq. in the presence 
of thermal fluctuations (VSz = Vg = 0). One might worry 
that this problem does not have a unique solution as both the 
left-hand and the right-hand side of Eq. 0 are not known 
a priori. Here one can, however, make use of the fact that 
Kubo’s fluctuation-dissipation theorerrP^ constraints the ther¬ 
mal forces on the skyrmion described by Fth in Eq. Q by 
linking them directly to the dissipative contributions of 
On average (Fth = 0), but its autocorrelation is proportional 
to the temperature and friction coefficients. In general it is 
given by 

(FJh(u;)F(Ja;')) = kBT[G^\u;) + GJ^ {-u;)]27rSico + J). 

(7) 

For small w one obtains (Fjjj(a;)F^|^(a;')) = 

AirksT aV 6{uj + oj') while off-diagonal correla¬ 
tions arise from the gyrodamping (Fj]^a;)F((t,(tLi')) = 
ATTiojksT aT 5{uj + uj'). Using Eq. 0 and demand¬ 
ing furthermore that the solution of Eq. 0 reproduces 
the correlation function {Ri(t)Rj{t')) (or, equivalently, 
{{Ri{t) — Rj{t'))'^)) obtained from the micromagnetic 
simulations, leads to the conditiorP^ 
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FIG. 2: Time dependence of the correlation function 

((Ri(fo-I-f) — Ri(fo))^) for T = O.IJ and different values 
of the Gilbert damping a (ujp = B — 0.0278 J is the frequency for 
cyclotron motion). 


on the (frequency-dependent) effective mass, gyrocoupling, 
damping and gyrodamping of the skyrmion. Furthermore, the 
corresponding spectrum of thermal fluctuations is given by 

Eq.Q__^ 

Fig.|^shows {{ARf)t = {{Rxih + t) - Rxito))"^). As 
expected, the motion of the skyrmion is diffusive; the mean 
squared displacement grows for long times linearly in time 
((Ai?)^)t = 2Dt, where D is the diffusion constant. Usu¬ 
ally the diffusion constant of a particle grows when the fric¬ 
tion is lowered^. For the skyrmion the situation is opposite: 
the diffusion constant becomes small for the small friction, 
i.e., small Gilbert damping a. This surprising observation has 
its origin in the gyrocoupling Q: in the absence of friction 
the skyrmion would be localized on a cyclotron orbit. From 
Eq. 0, we obtain 


D = ksT' 


oD 


+ {aVf 


(9) 


The diffusion is strongly suppressed by Q. As in most materi¬ 
als a is much smaller than unity while V ^ G, the skyrmion 
motion is characterized both by a small diffusion constant 
and a small friction. Such a suppressed dynamics has also 
been shown to be important for the dynamics of magnetic 
vorticeJ^. For typical parameters relevant for materials like 
MnSi we estimate that it takes 10“® s to 10“® s for a skyrmion 
to diffusive over an average length of one skyrmion diameter. 


To analyze the dynamics on shorter time scales we show in 
Fig .|^four real functions parametrizing G ^ (w): a frequency- 
dependent mass m{uj), gyrocoupling (/(w), gyrodamping 
ar(a;) and dissipation strength a'D{uj) with 


We therefore determine first in the presence of thermal fluc¬ 
tuations (VS^ = rig = 0) from simulations of the stochastic 
LEG equation 0 the correlation functions of the velocities 
and use those to determine Gijiw) using Eq. 0. After a sim¬ 
ple matrix inversion, this Axes the left-hand side of the equa¬ 
tion of motion, Eq. 0, and therefore contains all information 


„_i, s _ l'aT>{uj) — iujm{uj) —(/(w)-f zQ;a;r(u;)\ 

^ \ G{ij->) — iu}ar{uj) aD{u}) — iujm{u}) J 

For w —0 one obtains the parameters of Eq. 0. All pa¬ 
rameters depend only weakly on temperature, G and m are ap¬ 
proximately independent of a, while the friction coefficients 
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FIG. 3: Dissipative tensor aD, mass m, gyrocoupling Q and gyro¬ 
damping qF as functions of the frequency uj for the diffusive motion 
at T = 0.1 (solid lines). They differ strongly from the “apparent” 
dynamical coefficients (see text) obtained for the force driven (red 
dashed line) and current driven motion (green dot-dashed line). We 
use a = 0.2, /3 = 0.1. The error bars reflect estimates of systematic 
errors arising mainly from discretization effects, see appendix [B] 


bhT=0.05 wT=0.15 
AAT=0.1 <k>T=0.2 




FIG. 4: Dissipative strength aD, mass m, gyrocoupling Q and gy¬ 
rodamping aF as functions of the Gilbert damping a for different 
temperatures T. 


aV and ofF are linear in a, see Fig. In the limit T ^ 0, 
Q{uj ^ 0) tak es the value —Ttt, fixed by the topology of the 
skyrmiodSEi. 

Both the gyrodamping F and and the effective mass m 
have huge numerical values. A simple scaling analysis of the 
Landau-Lifshitz-Gilbert equation reveals that both the gyro¬ 
coupling Q and T> are independent of the size of the skyrmion, 
while F and m are proportional to the area of the skyrmion, 
and frequencies scale with the inverse area, see appendix 
[B] For the chosen parameters (the field dependence is dis¬ 


cussed in the appendix |^, we find m « 0.3 Nfnp uiq and 
aF « aO.TWflip mo, where mo = is the mass of a sin¬ 
gle flipped spin in a ferromagnet (1 in our units) and we have 
estimated the number of flipped spins, Nmp, from the total 
magnetization of the skyrmion relative to the ferromagnetic 
background. As expected the mass of skyrmions grows with 
the area (consistent with an estimatd^for m obtained from the 
magnon spectrum of skyrmion crystals), the observation that 
the damping rate aV is independent of the size of skyrmions 
is counter-intuitive. The reason is that larger skyrmions have 
smoother magnetic configurations, which give rise to less 
damping. For realistic system parameters J = 1 meV (which 
yields a paramagnetic transition temperature Tc 10 K, but 
there are also materials, i.e. FeGe, where the skyrmion lattice 
phase is stabilised near room-temperaturdi^ and a = 5 A and 
a skyrmion radius of 200 A one finds a typical mass scale of 
10-25 kg. 

The sign of the gyrodamping aF is opposite to that of the 
gyrocoupling Q. This implies that aF describes not damp¬ 
ing but rather anf/damping; there is less friction for cyclotron 
motion of the skyrmion than for the linear motion. The nu¬ 
merical value for the antidamping turns out to be so large 
that Vm + F^ < 0. This has the profound consequence that 
the simplified equation of motion shown in Eq. Q cannot be 
used: it would wrongly predict that some oscillations of the 
skyrmion are not damped, but grow exponentially in time due 
to the strong antidamping. This is, however, a pure artifact 
of ignoring the frequency dependence of G~^{uj), and such 
oscillations do not grow. 

Fig.j^shows that the dynamics of the skyrmion has a strong 
frequency dependence. We identify the origin of this fre¬ 
quency dependence with a coupling of the skyrmion coordi¬ 
nate to pairs of magnon excitations as discussed in Ref. [31] 
Magnon emission sets in for w > 2ijjp where ujp — B is the 
precession frequency of spins in the ferromagnet (in the pres¬ 
ence of a bound state with frequency uib, the onset frequency 
is u!p + ojb. Ref. EB. This new damping channel is most ef¬ 
ficient when the emitted spin waves have a wavelength of the 
order of the skyrmion radius. 

As a test for this mechanism, we have checked that only 
this high-frequency damping survives for a —0. In Fig. 
we show the frequency dependent damping a'D{uj) for various 
bare damping coefficients a. For small uj it is proportional to 
a as predicted by the Thiele equation. For uj > 2ujp, however, 
an extra dampling mechanism sets in: the skyrmion motion 
can be damped by the emission of pairs of spin waves. This 
mechanism is approximately independent of a and survives 
in the a —0 limit. This leads necessarily to a pronounced 
frequency dependence of the damping and therefore to the ef¬ 
fective mass m(u’) which is related by the Kramers-Kronig 
relation m(uj) = - aV(uj ) ^ aViuj). Note also that 
the large a independent mass m{uj 0) is directly related to 
the a independent damping mechanism for large oj. Also the 
frequency dependence of m{uj) and G{uj) can be traced back 
to the same mechanism as these quantities can be computed 
from q.'D{oj) and aF(w) using Kramers-Kronig relations. For 
large frequencies, the effective mass practically vanishes and 
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FIG. 6: Dynamical coupling coefficients for the force driven motion 
(a = 0.2). In the static limit everything but the real part of the diago¬ 
nal vanishes. Re Sg^ (to) however approaches the total magnetization 
Mtat as expected. The error bars reflect estimates of systematic er¬ 
rors, see appendix [B| 

the ‘gyrocoupling’ Q drops by a factor of a half. 


IV. FORCE-DRIVEN MOTION 

Next, we study the effects of an oscillating magnetic field 
gradient in the absence of thermal fluctuations. As 

the skyrmion has a large magnetic moment relative to the 
ferromagnetic background, the field gradient leads to a force 
acting on the skyrmion. In the static limit, the force is exactly 
given by 

Fg(cu ^ 0) = (10) 

Using G“^(w) determined above, we can calculate how the 
effective force Sg{u})'VB^iu}) (see Eq. depends on fre¬ 
quency. Fig.|^shows that for a; —0 one obtains the expected 


result Sg{uj —0) = while a strong frequency de¬ 

pendence sets in above the magnon gap, for uj > ojp. This 
is the precession frequency of spins in the external magnetic 
field. 

In general, both the screening of forces (parametrized 
by Sg(a;)) and the internal dynamics (described by 
G“^(a;)) determines the response of skyrmions, V{uj) = 
G{ cu) Sgi^uj^ V B^(^oj'j. Therefore it is in general not possi¬ 
ble to extract, e.g., the mass of the skyrmion as described by 
G~^{uj) from just a measurement of the response to field gra¬ 
dients. It is, however, instructive to ask what “apparent” mass 
one obtains, when the frequency dependence of Sg(a;) is ig¬ 
nored. We therefore define the “apparent” dynamics G~^{uj) 
by Ga{uj)Sg{uj = 0) = G(Lj)Sg(uj). The matrix elements 
of G~^(ui) are shown in Fig. ^as dashed lines. The appar¬ 
ent mass for gradient-driven motion, for example, turns out 
to be more than a factor of three smaller then the value ob¬ 
tained from the diffusive motion clearly showing the impor¬ 
tance of screening effects on external forces. The situation 
is even more dramatic when skyrmions are driven by electric 
currents. 


V. CURRENT-DRIVEN MOTION 

Currents affect the motion of spins both via adiabatic and 
non-adiabatic spin torque^. Therefore one obtains t wo ty pes 
of forces on the spin texture even in the static limiPSESlil 

The effect of a time-dependent, spin-polarized current on 
the magnetic texture can be modelled by supplementing the 
right hand side of eq. ([^ with a spin torque term Tst, 

Tst = -(v, ■ V)Mr + ^ [Mr X (v, ■ V)Mr] . (11) 

The first term is called the spin-transfer-torque term and is 
derived under the assumption of adiabaticity; the conduction- 
electrons adjust their spin orientation as they traverse the mag¬ 
netic sample such that it points parallel to the local magnetic 
moment Mr owing to Jb and Jsd- This assumptions is justi¬ 
fied as the skyrmions are rather large smooth objects (due to 
the weakness of spin-orbit coupling). The second so called /?- 
term describes the dissipative coupling between conduction- 
electrons and magnetic moments due to non-adiabatic effects. 
Both a and /3 are small dimensionless constants in typical ma¬ 
terials. From the Thiele approach one obtains the force 

Fg{uj —y 0) — ^ X -f I3T) Vg. (12) 

For a Galilei-invariant system one obtains a = /?. In this 
special limit, one can easily show that an exact solution of the 
LLG equations in the presence of a time-dependent current, 
described by Vs{t) is given by M{r — Vs(t')dt') pro¬ 
vided, M(f) is a static solution of the LLG equation for Vg = 
0. This implies that for a = /3, the skyrmion motion exactly 
follows the external current, R{t) = Vg{t). Using Eq. O, 
this implies that for a = /3 one has G~^{uj) = Sc{uj)- Defin¬ 
ing the apparent dynamics, as above, Ga{oj)Sc{uj = 0 ) = 
G{uj)Sciio) one obtains a frequency independent G~^{uj) = 
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FIG. 7: Dynamical coupling coefficients (symbols) for the current- 
driven motion (a = 0.2, /3 = 0.1, J = 1, A = 0.18 J, B = 0.0278). 
These curves follow almost the corresponding matrix elements of 
(w) shown as dashed lines. A deviation of symbols and dashed 
line is only sizable for Re 5^^. 


o —0 a=0.2, p=0 

□—□ a=0.2, p=0.1 


«— »a=0.2, p=0.15 
a=0.2, p=0.19 
>—> a=0.2, p=0.3 



FIG. 8: Mass m{uj) and gyrodamping ar(a;) as functions of the 
driving frequency ui for the current-driven motion. Note that both M 
and r vanish for a = (3. 


Sc{oJ = 0 ) = PVI — iUyQ: the apparent effective mass and 
gyrodamping are exactly zero in this limit and the skyrmion 
follows the current without any retardation. For a ^ f3, the 
LLG equations predict a finite apparent mass. Numerically, 
we find only very small apparent masses, oc a — /?, see 
dot-dashed line in upper-right panel of Fig. where the case 
a = 0.2, /3 = 0.1 is shown. This is anticipated from the anal¬ 
ysis of the a = P case: As the mass vanishes for a = P = Q, 
it will be small as long as both a and P are small. Indeed 
even for a 7 ^ /3 this relation holds approximately as shown 
in Fig. 1^ The only sizable deviation is observed for Re 
for which the Thiele equation predicts Re 0) = PV 

while Re —>■ 0) = oD as observed numerically. 

A better way to quantify that the skyrmion follows the cur¬ 


rent even for a 7 ^ /? almost instantaneously is to calculate 
the apparent mass and gyrodamping for current driven mo¬ 
tion, where only results for a = 0.2 and /3 = 0.1 have been 
shown. As these quantities vanish for a = P, one can ex¬ 
pect that they are proportional to a — /3 at least for small a, p. 
This is indeed approximately valid at least for small frequen¬ 
cies as can be seen from Fig. Interestingly, one can even 
obtain negative values for p > a (without violating causal¬ 
ity). Most importantly, despite the rather large values for a 
and P used in our analysis, the apparent effective mass and 
gyrodamping remain small compared to the large values ob¬ 
tained for force-driven motion or the intrinsic dynamics. This 
shows that retardation effects remain tiny when skyrmions are 
controlled by currents. 


VI. CONCLUSIONS 


In conclusion, we have shown that skyrmions in chiral mag¬ 
nets are characterised by a number of unique dynamical prop¬ 
erties which are not easily found in other systems. First, their 
damping is small despite the fact that skyrmions are large 
composite objects. Second, despite the small damping, the 
diffusion constant remains small. Third, despite a huge iner¬ 
tial mass, skyrmions react almost instantaneously to external 
currents. The combination of these three features can become 
the basis for a very precise control of skyrmions by time- 
dependent currents. 

Our analysis of the skyrmion motion is based on a two- 
dimensional model where only a single magnetic layer was 
considered. All qualitative results can, however, easily be 
generalized to a film with Nl layers. In this case, all terms 
in Eq. 0 get approximately multiplied by a factor Nl with 
the exception of the last term, the random force, which is en¬ 
hanced only by a factor y/Wp. As a consequence, the diffu¬ 
sive motion is further suppressed by a factor l/^/Wp while 
the current- and force-driven motion are approximately unaf¬ 
fected. 

An unexpected feature of the skyrmion motion is the an¬ 
tidumping arising from the gyrodamping. The presence of 
antidumping is closely related to another important property 
of the system: both the dynamics of the skyrmion and the ef¬ 
fective forces acting on the skyrmion are strongly frequency 
dependent. 

In general, in any device based on skyrmions a combination 
of effects will play a role. Thermal fluctuations are always 
present in room -temp erature devices, the shape of the device 
will exert force^SEl and, finally, we have identified the cur¬ 
rent as the ideal driving mechanism. In the linear regime, the 
corresponding forces are additive. The study of non-linear 
effects and the interaction of several skyrmions will be impor¬ 
tant for the design of logical elements based on skyrmions and 
this is left for future works. As in our study, we expect that 
dynamical screening will be important in this regime. 
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One therefore needs a method which focuses only on the 
region close to the skyrmion center. To locate the skyrmion, 
we use the z-component of the magnetization but take into ac¬ 
count only points where Ms(r) < —0.7 (the magnetization of 
the ferromagnetic background at T = 0 is -t-1). We therefore 
use 


p(r) = (1-M,(r))0[-M,(r)-O.7] (A3) 

where 0[a;] is the theta function. A first estimate, Rest = Rv, 
for the radius is obtained from 


FIG. 9: Skyrmion density based on the normalized «-component of 
the magnetization. 


J^rpir)d^r 

^ Ia 


(A4) 
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by integrating over the full sample volume V. Rest is noisy 
due to the problems mentioned above but for the system 
sizes simulated one nevertheless obtains a good first esti¬ 
mate for the skyrmion position. This estimate is refined by 
using in a second step for the integration area only D = 
{r € I |r — Rost I < r} where r is choosen to be larger 
than the radius of the skyrmion core (we use r = 1.3y/ 
where is the number of spins with < —0.7). Thus 
we obtain a reliable estimate, R — Rd, not affected by spin 
fluctuations far away from the skyrmion. Prom the resulting 
R{t), one can obtain the velocity-velocity correlation function 

{V,{to+t)Yj{to)). 

Appendix B: Scaling invariance 


Appendix A: Definition of the Skyrmion’s centre coordinate 

In order to calculate the Green’s function, Eq. (3), one 
needs to calculate the velocity-velocity correlation function. 
Therefore it is necessary to track the skyrmion position 
throughout the simulation. Mostly two methods have been 
used so far for thiJ^ (i) tracking the centre of the topological 
charge and (ii) tracking the core of the Skyrmion (reversal of 
magnetization). 

The topological charge density 

Ptop(i’) = ;pn(r) • (5::cn(r) x dyh{r)) (Al) 

47r 

integrates to the number of Skyrmions in the system. There¬ 
fore for our case of a single Skyrmion in the ferromagnetic 
background this quantity is normalized to 1. The center of 
topological charge can therefore be defined as 

R = J (frptop{r)r (A2) 

Por the case of finite temperature this method can, however, 
not be used directly. Thermal fluctuations in the ferromagnetic 
background far away from the skyrmion lead to a large noise 
to this quantity which diverges in the thermodynamic limit. 
A similar problem arises when tracking the center using the 
magnetization of the skyrmion. 


To obtain analytic insight into the question of how param¬ 
eters depend on system size and to check the numerics for 
lattice artifacts it is useful to perform a scaling analysis of the 
sLLG equations and the effective equations of motion for the 
skyrmion. 

We investigate a scaling transformation, where the radius 
of the skyrmion is enlarged by a factor p, M(r) -> M(r) = 
M(r/p). 

Por the scaling analysis, we use the continuum limit of 
Eq. (4) (setting a = 1) 

^ (VM)^ -l-AM-VxM-B-M 

The three terms scale with 77 °, rj and respectively. To ob¬ 
tain a larger skyrmion, we therefore have to rescale A —> A/p 
and B —> B/p^. This implies that the Bed term in the sLLG 
equation scales with 1 /p^ and therefore also the time axis has 
to be rescaled, t -A rft, implying that all time scales are a 
factor of p^ longer and all frequencies a factor 1 /rf smaller. 
Similiarly, the driving current, i.e. Vg is reduced by a factor 
p“^ while the temperature remains unsealed. This implies that 
when M(r, t) is a solution for a given value of A, B and Vg 
and G( oj) the corresponding velocity-correlation function of 
the skyrmion, then M(r/p, t/rf) is a solution for A/p, Bjrf', 
Vgjr] with correlation function G(a;p^). 

Accordingly, the w —0 limit and therefore the gyrocou¬ 
pling Q, the friction constant aD and the diffusion constant of 
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FIG. 10: Field dependence of the dissipative tensor aT>, the mass m, 
the gyrocoupling Q and the gyrodamping aF for J = 1, A = 0.18 J. 
The main effect is that both m and qF shrink when the size of the 
skyrmion shrinks with increasing ■ 

the skyrmion are independent of p, consistent with the analyt¬ 
ical formulas eq. (6). In contrast, the mass of the skyrmion, 
TO, and the gyrodamping aF scale with They are therefore 
proportional to the number of spins constituting the skyrmion 
consistent with our numerical findings. When one does, how¬ 
ever, change the external magnetic field fox fixed strength A 
of the DM interaction, the internal structure of the skyrmion 
and therefore also G(w) change quantitatively in a way not 
predictable by scaling. In Fig. [^we therefore show the in¬ 


dependence of aV, TO, Q and aF. For increasing B the size of 
the skyrmion shrinks. From our previous analysis, it is there¬ 
fore not surprising that both the effective mass to and the gy¬ 
rodamping aF shrink substantially while oD and Q are less 
affected . Quantitatively, this change of to and aF is how¬ 
ever not proportional to the number of spins constituting the 
skyrmion. 

We have tested numerically the scaling properties for rj = 
1.2 and find that all features are quantitatively reproduced. 
Small variations on the level of a few percent do, however, 
occur reflecting the typical size of features arising from the 
discretization of the continuum theory. A conservative esti¬ 
mate of such systematic discretization effects for the diffusive 
motion is given by the error bars in Figs. (all statistical er¬ 
rors are smaller than the thickness of the line). For the field- 
driven motion (Fig. |^and Fig. spatial discretization effects 
lead to a different source of errors. For very small held gradi¬ 
ents and high frequencies the displacement of the skyrmion is 
much smaller than the lattice spacing and the response is af¬ 
fected by a tiny pinning of the skyrmion to the discreet lattice. 
For larger gradients, however, nonlinear effects set in and for 
small frequencies the skyrmion starts to approach the edge of 
the simulated area. In Fig. we therefore used for the force- 
driven motion \7B = 0.0005 for w < 2(jjp and VB = 0.0015 
for oj > 2ujp. Error bars have been estimated from variations 
of the numerical values when VB was varied from 0.0001 to 
0.0015. For the current-driven motion errors are so tiny that 
they are not shown. 
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